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Q\ , Abstract: The single rung t-J ladder is analyzed in a mean field theory using 



Gutzwiller renormalization of the matrix elements to account for strong correlation. 
2 ■ The spin liquid (RVB) state at half-filling evolves into a superconducting state upon 

'-^ ■ doping. The order parameter has a modified d-wave character. A lattice of weakly 

coupled ladders should show a superconducting phase transition. 
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There is a striking difference between the properties of a chain and a ladder (dou- 
ble chain) antiferromagnetic (AF) s=l/2 Heisenberg model. Whereas the chain has 
power law decay of the AF-correlations, the ladder has a purely exponential de- 
cay and a finite energy gap in the spin excitation spectrum, i.e. a spin gap (see 
for example Ref.l). If, as is the case for other spin gap systems, the spin gap 
persists to finite doping, then the possibilities for superconducting fluctuations are 
greatly enhanced in ladder systems [2] . Recently we pointed out that the compound 
Sr2Cu406 offers the possibility of realizing a lattice of weakly coupled ladders [3,4]. 



This compound is a member of the homologous series Sr„_iCu„_|_i02„, which differ 
from known high-Tc cuprates through the presence of a parallel array of line defects 
in the Cu02-planes. In this letter we examine the properties of undoped and lightly 
doped ladders described by a t-J model within a mean field approximation which 
uses a Gutzwiller renormalization factor to approximate the local constraint. Super- 
conducting correlations of a modified d-wave symmetry are predicted. The system 
of weakly coupled chains offers an interesting example of a short range resonance 
valence bond (RVB) state in a system intermediate between 1 and 2 dimensions [5] . 
The basic model we apply is the t-J model on the ladder with the Hamiltonian 

ri = — 2_^V'x /_^ '^i+l,a,s'^i,a,s + ty'^i,l,s'^i,2,s + tl.C.\ 
i,s a=l,2 

+ 2^1'^^ A^ ^i+l,a ■ ^i,a + Jy'^i,! ' ^j,2j 
i a=l,2 

where i runs over all rungs and a over the two legs, 1 and 2. The constraint 
Ss Ci,a,sCj,a,s < 1, projcctiug out all the doubly occupied states on each site, is 
implied. We examine this model within the mean field theory used by Zhang et 
al. for the 2D square lattice [6]. The constraint is taken into account approx- 
imately by a Gutzwiller type renormalization of the coupling constant along the 

Cham, \Cj_|_;^ (j gCj^a^s) = 9tx\Ci+l^a,s'^i,a,s/0i ^-^d (Oj_|_i^(j • Oi^a) = 9jx\'^i+l,a ' '^i,a/Oj ^^'i 

analogous gty and gjy on the rungs with (...) and (...)o the expectation values in the 
projected and unprojected state respectively. These (yf-factors are determined by the 
ratio of the probabilities of the matrix elements in the projected and unprojected 
states [6]. The renormalized coupling constants are then defined by Ja = JaQja and 

Using the fermion representation for the Heisenberg term in (1), we introduce 
two types of mean fields, Xx = (cl+i,a,sCi,a,s)o and A^ = {ci+i,a,iCi,a^)o, and Xy = 
(cIi,sCi,2,s)o and Ay = (cj_i^j^Cj_2t)o) a bond and pairing mean field, respectively. The 
t-J model is then reduced to a renormalized Hamiltonian with the constraint released 



T-Cmf = X!^kCk,sCk,s - X![^kCkiC_kT + h.c] 

k,s k 

(2) 

+^ E /.^«(ix«r+iA.r) 

Q!=a;,j/ 

where /a; = 2 and fy = l and A^ is the number of lattice sites. The summation over 
k is restricted to — tt < k,j. < -k and ky equals to (bonding) or tt (antibonding). In 

(2) 

3 ~ 

£k = -2 E /"('^o + T'^aXa)cOsA;Q - yU 
a ^ 

3^ ~ (3) 

Ak = T^X. faJa^aCOska 
^ a 

with /i as the effective chemical potential. We calculate the Gutzwiller renormaliza- 
tion factors by including the correlations of the probability between the two nearest 
neighbor sites [7], which improves the approximation of Ref.6. They are given by 

g^^ =2n6l{n{l + 6)+Axl) 

gj^ = 4nV(n^(l + ^f + ^5\\K? - \Xa?) + 16(|A«r + \x.\')) ^^^ 



where n is the electron density and 5 = 1 — n is the hole doping concentration. 

The mean field Hamiltonian has to be solved self-consistently together with the 
mean field dependent Gutzwiller renormalization factors. Additionally, the chemical 
potential is chosen to give the correct electron concentration, n = 1 — 6 = (nj)o. 
Numerical solutions of these mean fields as a function of doping concentration are 
plotted in Fig.l for realistic parameters, tx = ty and Jx = Jy = O.Stx- Within 
this mean field treatment the excitations are described by the effective quasiparticle 
Hamiltonian obtained by standard Bogolyubov transformation of Ti-MF, T^eff = 
J2k,a=i,2 -E'k7ko-7ko- with the spectrum E'^ = 6"^ + A^. The quasiparticles define the 
the ground state by the condition 7k,o-|'^o) = 0. For finite Ak these excitations have 
always a finite gap. 



We begin at half-filling where this system reduces to a Heisenberg ladder. This 
case was considered by various groups for different couplings along the chains, J^, 
and rungs Jy [1,8,9]. In the hmit Jy ^ Jx > the ground state consists essentially 
of singlet dimers on each rung leading to a spin liquid (singlet) ground state. The 
lowest spin (triplet) excitation is obtained by replacing one singlet by a triplet dimer 
with an excitation energy ~ Jy — Jx- By numerical diagonalization and dimer mean 
field treatments it was demonstrated that the spin liquid ground state persists even 
for Jx ~ Jy with a spin gap of 0.4 — 0.6 J for J = Jx = Jy [8,9]. 

The mean field solution at half-filling independent of the coupling constants 
satisfies the relation A^-Aj^ + XxXy = 0. This state is similar to the d-wave RVB 
in the square lattice, where the pairing mean fields A^ differ by a phase tt in x- 
and ^/-direction. It is also identical to the Affleck-Marston flux phase [10] with half 
integer flux quanta as follows from a SU(2) symmetry [6]. Upon doping, the flux 
phase and the d-wave RVB state differ. In the case of the square lattice, the d-wave 
RVB state is found to have lower energy than the staggered flux phase [11]. We 
expect the similar result for the ladder, and will not discuss the flux phase here. 

Turning to the spin excitation (singlet-triplet) we find that it is not well described 
by the simple quasiparticle Hamiltonian Tie// given above. The spin gap is generally 
too large (^ 2 J for Jx = Jy). Additionally the excitation energy has several minima 
at the momenta q = (0, 0), (0, tt), (tt, 0) and (tt, tt), while in the theories mentioned 
above the only minimum lies at q = (tt, tt) [8,9]. 

We can remedy this flaw in our treatment partially by taking the residual inter- 
action among quasiparticles into account. Most easily the spectrum of the spin 
excitations can be obtained from the dynamical (transverse) spin susceptibility, 
x(u;,q) = {{SZf^, S^))^. We use the equation of motion to determine x(c<j,q) on 
the level of RPA including the residual interaction. In Fig.2 we show the spectrum 
Imxi^, q) at q = (vr, tt) for both the simple mean field and the RPA corrected result 
with Jx = Jy The RPA shows a sharp excitation peak at ^ 0.15 Jx and a broad 
continuous spectrum between about 2Jx - 3Jx which is the strongly suppressed re- 
mainder of the continuous excitation spectrum of Tie//- The analysis of the full 
q-dependence shows that the minimal excitation is obtained at q = (7r,7r). The 



corresponding binding energy is rather large bringing the excitation gap down to a 
value of 0.15 Jj: which is only slightly smaller than the one found in other calculations 
(^ 0.5J,) [8,9]. 

As a consequence of the local constraint the singlet dimer on a rung is a coherent 
superposition of pair states in the bonding (q_|_ s = (q i s + Ci,2,s)/v^) ^^'^ the anti- 
bonding (q_ s = (4,1,8 ~ w,2,s)/v^) state, i.e.: avoiding double occupancy the rung 
state (l/v^)(cj^_|ct^| — c|„|c|_|)|0) gains the maximal exchange energy. In other 
words, the splitting of the bonding and antibonding band, which would appear in 
a band structure description, is absent in the present of the local constraint so that 
they are equally filled at half-filling. 

To get an intuitive understanding of the properties of this sharp excitation found 
in RPA let us consider the problem in the strong coupling limit, Jy ^ J^- In the 
mean field a triplet excitation creates two quasiparticles which may be considered 
as spinous freely moving and destroying each one singlet dimer. The energy loss is 
3Jy/2, neglecting the kinetic energy contribution. These two quasiparticles can form 
a triplet dimer on a rung with the lower energy Jy. Thus, there exists an effective 
attraction between them leading to a bound exciton state within the quasiparticle 
excitation gap. This interaction has only an attractive triplet channel if one of 
the two quasiparticles is in the bonding and the other in the antibonding band. 
Furthermore, the kinetic energy of the triplet dimer is given by JyCosqx which is 
lowest at Qx = TT. Hence, we interpret the sharp excitation in the RPA calculation 
for Jx = Jy as an exciton with lowest energy at q = (tt, tt) and the reduced continuum 
above ~ 2Jx as a spinon continuum (Fig.2). 

Upon doping holes the two quasiparticle bands which are degenerate at half-filling 
split. The antibonding band is raised and the bonding band slightly lowered relative 
to the chemical potential. As holes are doped the kinetic energy can be lowered if 
they occupy preferentially the antibonding band. In this way the antibonding band 
is gradually emptied down to a critical doping 6c where eventually only the bonding 
band is occupied. The loss of antibonding electrons leads to a gradual decrease of 
the pairing amplitude Aq, which disappears eventually at Sc- The long tail of finite 
Aq, down to 6c may be an artifact of our approximation. 



The behavior of A^^ yields a monotonic decrease of the quasiparticle excitation 
gap described by Tie//- On the other hand, the spin excitation band formed by 
the triplet exciton is shifting non-monotonically with doping leading to an increase 
in the spin gap at small doping (see Fig.3, for coupling constants t^ = ty = t and 
Jx = Jy = J = 0.3t). With the shift of the Fermi levels in the antibonding band the 
optimal relative momentum q^ for a triplet exciton of two quasiparticles deviates 
gradually from vr. Since Qx = t^ yields the strongest attraction this deviation weakens 
the interaction and decreases the binding energy. The resulting increase of the 
spin gap for small doping shows that here hole doping suppresses antiferromagnetic 
correlation and stabilizes a spin liquid state. Also results of exact diagonalization 
confirm this tendency [12]. 

The BCS superconductivity order parameter is given by the pairing mean field 
multiplied by the Gutzwiller renormalization factor gta as shown in Ref.6. In Fig. 4 
we show the two components of Asca = ^aQta- There is a phase difference of vr be- 
tween Ascx and A^^^, reminding of d-wave pairing, although d-wave is a misnomer 
in this low symmetry system. 

In summary, our mean field theory of the t — J-ladder model gives a description of 
the doped spin liquid system. The spin liquid state persists for a finite doping region 
away from half-filling and seems even to become more stable with weak doping. 
Although mean field is certainly not a good concept for quasi ID system, we have 
seen that qualitatively reasonable results are obtained. We may expect that weak 
interchain (intraplanar and interplanar) coupling would lead essentially to a 3D 
situation stabilizing the mean field solution. The superconductivity obtained in this 
theory is intimately connected with the existence of a spin liquid state (see also in 
Ref.2 and 13). Considering their energy scales we observe that at light doping the 
spin gap is by far larger than the superconducting pair correlation energy. Beside 
the triplet exciton mode discussed here we expect in the doped region a low lying 
collective (sound) mode in connection with the superconducting order will occur 
[14]. 

Our result suggests that also for chains extended in ^/-direction having 4, 6 or a 
larger even number of parallel chains would qualitatively be very similar. Although 



the spin gap must be diminished with growing system size, as recent mean field 
calculations show [9], the concept of the doped spin liquid system is still applicable. 
In this sense the superconducting state would continuously tend to a real d-wave 
state, if we would extend the lattice in y-direction approaching in this way the 2D 
square lattice. 
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Figures Caption 



Fig.l: The mean fields as a function of tlie doping concentration for J^ = Jy = J 
0.3t. 



Fig. 2: Tlie spectrum of tlie spin-excitations for J^ = Jy = J = 0.3t and q = (tt, tt). 
Tlie solid line denotes the exciton spectrum obtained from RPA and the dashed 
line the mean field quasiparticle spectrum neglecting the interaction among the 
quasiparticles. The sharp exciton mode at cj ^ 0.15 J contains practically all weight 
of the excitation spectrum. 



Fig. 3: The spin gap as a function of the doping concentration for J^^ = Jy = J = 0.3t. 
The solid line denotes the bottom of the quasiparticle spectrum and the dashed line 
the gap of the exciton state. 

Fig. 4: The superconductivity order parameter versus the doping concentration for 



